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Analytic solution of the fractional advection diffusion equation for the time-of-flight 

experiment in a finite geometry 

B.W. PhilippaQ R-D. White, and R.E. Robson 
School of Engineering and Physical Sciences, James Cook University, Townsville 4811, Australia 

A general analytic solution to the fractional advection diffusion equation is obtained in plane 
parallel geometry. The result is an infinite series of spatial Fourier modes which decay according 
to the Mittag-LefBer function, which is cast into a simple closed form expression in Laplace space 
using the Poisson summation theorem. An analytic expression for the current measured in a time- 
of-flight experiment is derived, and the sum of the slopes of the two respective time regimes on 
logarithmic axes is demonstrated to be —2, in agreement with the well known result for a continuous 
time random walk model. The sensitivity of current and particle number density to variation of 
experimentally controlled parameters is investigated in general, and the results applied to analyze 
selected experimental data. 

PACS numbers: 05.40.Fb, 73.50.-h, 05.60.-k 
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I. INTRODUCTION 



Modern solid state electronics is largely based upon in- 
organic, crystalline materials, such as silicon and germa- 
nium, the transport properties of which are generally well 
understood The same applies to gaseous electronics, 
for which there is a one-to-one correspondence with crys- 
talline condensed matter 0. On the other hand, organic 
semiconductors are attracting increasing interesting be- 
cause of their desirable properties, such as transparency, 
flexibility, and the prospect of economic advantage over 
inorganic electronics [3|. Organic materials, which may 
be amorphous, exhibit electrical properties which are 
generally qualitatively and quantitatively quite different 
from inorganic materials [4]. For example, charge carriers 
in a time-of-flight experiment exhibit long lived, spatially 
dispersed structures. Furthermore, the roles of the mo- 
bility and diffusion coeflScients, ^ and respectively, 
are not at all clear cut, as they are in crystalline struc- 
tures or gases. Such anomalous or "dispersive" behavior 
arises because the scattering of charge carriers may be 
accompanied by trapping in localized states for times r, 
as determined by a "relaxation function" 0(t), which has 
an asymptotic time dependence ~ t~'^, with fractional 
exponent 7. 

The recent interest in "fractional kinetics" derives 
mainly from the seminal paper of Scher and MontroU 5] , 
whose discussion in terms of a continuous time random 
walk has spawned an extensive literature in its own right 
01 l6l-[ll| . In this literature, it is often assumed that the 
charge carrier number density n(z,t) may be found as 
the solution of a fractional diffusion equation, which for 
present purposes we will refer to as the "Caputo" form 



of the fractional advection diffusion equation: 
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where q Dj is the Caputo fractional partial derivative 
with respect to t of order 7. The Caputo derivative (see 
Appendix A) accounts for trapping in localized states. 
This is appropriate for a thin sample of amorphous ma- 
terial conflned between two large plane parallel bound- 
aries, with all spatial variation confined to the normal 
direction, which defines the z axis of a system of coor- 
dinates. In addition it is assumed that the small signal 
limit prevails, and that both the drift velocity W — jiE 
, also directed along the z axis, and the longitudinal dif- 
fusion coefficient derive entirely from an externally 
applied field E. For non-dispersive transport, 7—^1 
such that ^Djn ^ ^, and Eq. ([J) assumes the familiar 
classical form [l^j. The present article focuses on new 
techniques for solution of Eq. ([TJ for the purposes of 
better understanding the factors influencing experiment. 

Before proceeding with the detailed analysis, it is im- 
portant to bear in mind that Eq. ([T]) is only approximate. 
Just as the kinetic theory of classical charge carrier trans- 
port in crystalline semiconductors and gases has been de- 
veloped to a sophisticated level through solution of Boltz- 
mann's kinetic equation, a more general and accurate pic- 
ture of anomalous transport in amorphous media should 
be obtained through solution of a fractional kinetic equa- 
tion in phase space, in which the microscopic collision 
operator accounts for scattering and trapping processes. 
Projection onto conflguration space is achieved by in- 
tegration over velocity space, yielding (with approxima- 
tions) Eq. ([!]) plus expressions for macroscopic properties 
such as /i and D^. The phase space approach is beyond 
the scope of the present work, and the reader is referred 
to [13] for such considerations. 

Whatever the medium, gaseous or condensed matter, 
crystalline or amorphous, the advection diffusion equa- 
tion (dl is usually assumed to provide the link between 
theory and experiment, its limitations not withstanding. 
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Thus, on the one hand, solution of the Boltzmann ki- 
netic equation provides theoretical values of /i and 13^, 
and on the other, solution of Eq. ([T|) for n(z, t), with ap- 
propriate boundary and initial conditions, enables exper- 
imental data to be unfolded to furnish empirical values 
of the same transport properties. Comparison of theo- 
retically derived and experimentally measured transport 
properties then gives information about the fundamen- 
tal microscopic nature of the interaction of charge carri- 
ers with the medium including the trapping/detrapping 
process. This procedure is standard for electrons and 
ion "swarms" in gases [l2|, but application of the idea 
to amorphous media awaits the further development of 
fractional Boltzmann phase space kinetics. That is part 
of our long term theoretical program, but in the mean- 
time, we focus in the present article on the more practical 
imperative of developing an accurate and efficient means 
of solving Eq. (P). 

To this end, a simple and numerically efficient solution 
of Eq. dJ) would be highly desirable. Previously reported 
solutions of fractional diffusive systems for bounded me- 
dia have been expressed in terms of infinite series solu- 
tions [1, [lil, |T^. We show that the series solution to 
Eq. (P) with absorbing boundaries may be collapsed 
into a simple closed form solution in Laplace space by 
building upon the experience gained in solution of the 
non-dispersive diffusion equation in gaseous electronics, 
specifically, for the pulsed radiolysis drift tube experi- 
ment 0. The structure of this article is as follows: In 
Section II, we model the time of fiight experiment [l^ 
and obtain a formal analytic solution of Eq. ([T]) as a series 
of Mittag-Leffler functions, which is cast into a tractable 
form, suitable for practical purposes, using the Poisson 
summation theorem. In Section III, we express the cur- 
rent measured in a time-of-flight experiment in terms of 
this analytic solution, and show analytically that sums 
of the slopes in distinct time regimes add up to -2 on a 
log-log plot, as first predicted by Scher and Montroll [1] 
and as observed in many experiments In Section IV, 
we explore the way that current varies with experimental 
parameters, and go on to fit selected experimental data. 
We show that our solution demonstrates the power-law 
decay characteristic of dispersive transport. 



II. ANALYTIC SOLUTIONS OF THE 
FRACTIONAL DIFFUSION EQUATION 

In this article, we will use Eq. ([T]) to model a disor- 
dered semiconductor in a time of fiight experiment . 
The relationship between the various forms of the frac- 
tional advection diffusion equation using both Caputo 
and Riemann-Liouville forms of the fractional derivative 
operator are discussed in Appendix A. A one dimensional 
equation, such as ([1]), is appropriate for a thin sample 
of disordered material confined between two large plane 
parallel boundaries, which we shall take to be at z = 
and L respectively. All spatial variation is confined to 



the normal direction, which defines the z axis of a sys- 
tem of coordinates. In addition it is assumed that the 
small signal limit prevails, and that both the drift ve- 
locity W = iiE (where /i is the mobility) and the lon- 
gitudinal diffusion coefficient derive entirely from an 
externally applied field E. 

In the idealized time-of-flight experiment, a sharp pulse 
of no charge carriers is released from a source plane z = zq 
at time t = to, i.e.. 



n{z,to) = noS{z - zq). 



(2) 



and the fractional advection diffusion equation is solved 
using the methods and techniques described below. The 
solution for other experimental arrangements, e.g., for 
sources distributed in space and/or emitting for finite 
times, can be found by appropriate integration of this 
fundamental solution over zq and/or respectively. The 
solution for perfectly absorbing boundaries, for which 



n(0,t) = = n{L,t) 



(3) 
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where the spatial modes are 



it-tor), (4) 



(cos [km{z - Zq)] - COS [/c„(z + Zq)]) , 



and where 



= Dl (A^ + kl, 
rmr 



(5a) 
(5b) 
(5c) 



In Eq. (Ill), E-y{z) is the Mittag-Leffier function of order 
1- 



fc=0 

E^{z) = Ea.iiz) 



T{ak + f3) 
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Equation Q gives an exact solution, however, this ex- 
pression is somewhat difficult to manipulate due to the 
presence of the Mittag-Leffler function. Furthermore, a 
large number of terms are needed for this series to con- 
verge, and the numerical evaluation of the Mittag-Leffler 
to suitable precision is computationally difiicult. 

As is well known, fractional models obey a correspon- 
dence principle, where non-fractional behavior is recov- 
ered in appropriate limits. In this case, in the limit 7 — >■ 1 
the Mittag-Leffier function reduces to an exponential, i.e. 
Ei{z) = e^ and (gD reduces to Ec^ (3b) in Ref. ^. In 
the classical, non- fractional limit f2!|, it was shown that 
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the series convergence could be substantially improved 
through application of the Poisson summation theorem 
(PST): 



E 
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where F{k) is the Fourier transform of f{x). This arti- 
cle will demonstrate that the PST can also be applied 
to the fractional advection diffusion equation with sim- 
ilar benefits. Attempting to apply the PST directly to 
Eq. ^ results in an intractable Fourier transform in- 
volving the Mittag-LefHer function. On the other hand, 
the Mittag-LefHer function has a simple Laplace domain 
representation. Transformed into Laplace space, Eq. ^ 
becomes 



n(z,s) = no 



7-1 



m — 1 



(8) 



where without loss of generality we have taken to = 0. 

Applying the Poisson summation theorem to Eq. (jl 
gives the equivalent form 



n(z, s) 



ae 
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(9) 



where the space-independent parameters a and j3 are de- 
fined as 



a{s) 



(10) 



(11) 



Simpliiying Eq. (|9]), we obtain the closed form expres- 
sion 



n{z, s) 



^-/3|z-2n| 



4 sinh {f3z) sinh {(3zq) 



(12) 



A necessary condition for convergence to the closed form 
expression Eq. ([T^ is 



|exp(-2/3L)| < 1, 



(13) 



which defines the region of convergence of the Laplace 
domain function Eq. (jl2p . 

It should be emphasized that Eq. ([T2|) is a gen- 
eral result, valid for fractional and non-fractional cases. 
For normal transport (i.e., crystalline semiconductors or 
gaseous electronics), 7 = 1, and Eq. ([9]) has an ana- 
lytic inverse Laplace transform that reduces to Eq. (7) 
of 0, where it was obtained using time domain meth- 
ods. For dispersive transport, 7 < 1, and an analytical 
inverse Laplace transform is difficult to find, so the ap- 
plications presented below required numerical inversion 
of the Laplace transform[25|. 



III. CURRENTS AND THE SUM RULE 

A. Number, number density and charge carrier 
current in the time of flight experiment 

A typical time of fiight experiment measures the exter- 
nal current as photogenerated carriers are driven through 
the sample by an applied electric field. Under the condi- 
tion that the experimental time scale is much less than 
the RC time of the measurement circuit, the observed 
current is the space averaged conduction current 



(14) 



Expressed in terms of the number density n(z,i), the 
photocurrent is 



n{z,t)dz}, (15) 



where q is the charge on each carrier. The origin of Eq. 
(ITSI) is detailed in Appendix |B] Substituting the time 
domain n{z, t) solution Eq. Q into Eq. ((T5|) . the current 
is found to be 



I{t) = ^ Krat ^£^7,0 {-Ul^f ) 



(16) 



with 



2qnoe-^''>kmDL . , 

Km = ^^7-^ Sm [krnZo) 



[2XDl {e^"^ i-ir - 1) - Lujm] 



Alternatively, a closed form expression may be found 
in Laplace space by substituting Eq. [12] into Eq. (fT5|) . 



B. Sum rule for asymptotic slopes 

Experimental time of fiight current traces plotted on 
double logarithmic axes often demonstrate two distinct 
straight line regimes (see, for example. Figure [5]), a dis- 
tinctive shape which has been described as the "signa- 
ture" of dispersive transport In many materials, the 
sum of the slopes on logarithmic axes of these two regimes 
is very close to —2 (Refs. 0,[l3l), a prediction originally 
made for a continuous time random walk model by Scher 
and Montroll Q. In what follows, we prove that our ex- 
pression for the current, Eq. (1161) . demonstrates the same 
"sum of slopes" criterion. 

The small argument asymptote of the Mittag-Leffier 
function can be written down from its power series defi- 
nition, Eq. The result is 

£7,0 (-Wmi'^) -UJmt'', 



4 



where we have neglected terms of order O (^[ujm r]^) and 

higher. Substituting this into Eq. ([T6l) we find the early 
time current to be 

oo 

-^carly(^) ~ ^ ^ ^ra^ ^m^^ ^ 
m— 1 

Conversely, for the long time current, we use the large \z\ 
asymptote valid for negative real z [1^ 



z-^ 



+ 0[\z 



i-i-p 



If t is large, then by taking p = 1 we obtain the following 
form for the long time current 



_1 { -UJ,nt^) ^ 

r(-7) 



In summary, the asymptotic forms of the current for 7 ^ 
1 are 

ft-(i-T), early times 
'^^^^1t-(i+7), late times, ^ ^ 



in agreement with the sums of slopes condition. 

It is noteworthy that these asymptotes are indepen- 
dent of the boundary conditions imposed on the system. 
When solving the fractional diffusion equation n(z,t) is 
assumed to be factorable as n{z,t) — Z{z)T{t). The 
time-dependent function, T{t) can be expressed in terms 
of Mittag-Leffler functions 



t) 



(18) 



where Cm are the separation eigenvalues found by apply- 
ing the boundary conditions to the differential equation 
for Z{z). The asymptotes of the Mittag-Leffler functions 
[isj are such that physically acceptable solutions must 
have Cm < so that n{z,t) remains bounded as t — >■ 00. 
Imposing only the requirement that the boundary condi- 
tions result in a negative separation constant, using Eq. 
(|15p the current must take the form 



I{t) 



Using the asymptotic limits detailed above, the time de- 
pendence may be brought outside the summation, and 
the the same temporal asymptotes detailed above then 
follow. This result is independent of the spatial boundary 
conditions and hence independent of the specific form of 
Z{z). 



C. Transit Time 

The transit time can be obtained from the expres- 
sion for the total number of charge carriers within the 
medium. Defining 



Nis) 



we find in Laplace space 



n(z, s)dz 



N 



no 



-(A+/3)zo _ sinh(/3zo) , _ 

sinh {(3L) ^ 



(19) 

To simplify the mathematics and obtain an estimate for 
the transit time, we neglect diffusion by taking the limit 



NDr^=Q = ^ ( 1 - exp 



^s-' (L - zo) 



W 



(20) 



In the classical case with 7=1, the above equation 
has the expected inverse Laplace transform 



N 



(classic 
Dl=0 



[t) = Uq 



1- H t- 



W 



where H{t) is the Heaviside step function. 

For the dispersive case, where 7 < 1, Laplace inversion 
by complex contour integration gives 



NDL=o{t) = no 77, 



(21) 



where 



_ (—1)'"''""'^ sin (m7r7) F (7771) 



In the special case of 7 = 1/2, the power series Eq. ([2T 
is equivalent to the closed form expression 



<=°o''W-«oerf 



L - zo \ 
2WVi) 



(22) 



where erf is the Gaussian error function. It is interesting 
to note that Eq. ([2^ demonstrates great dispersion de- 
spite it being a zero diffusion limit of the true behavior 
of the system. 

A clear transit time cannot be precisely defined be- 
cause the packet of charge carriers becomes widely dis- 
persed. Nevertheless, there exist two regimes of cur- 
rent transport behavior, and the boundary between these 
regimes defines a "transit time" for the material. It can 
be seen that two distinct regimes will emerge from Eq. 
(mi, according to the magnitude of the term in paren- 
thesis. The transit time, defining the transition between 
regimes, is therefore approximately given by 



Zo 



1. 



5 



Solving for the transit time tf 



W 



1/7 



(23) 



This is in agreement with the expected experimental 
length and field dependence [1, [l3| ■ 
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Figure 1: (Color online) Impact of the fractional order 
7 on the temporal current profiles. Each curve is the 
current resulting from the respective number density 
solution of Figure [21 



IV. RESULTS 

A. Impact of model parameters on the density and 
current profiles 

The model discussed above has five parameters: the 
fractional drift velocity W, the fractional diffusion coef- 
ficient Dl, the fractional order 7, the initial source loca- 
tion zo, and the length of the sample L. These param- 
eters are constrained such that 0<7<l,0<zo<L 
and Dl > 0. The effects of varying the first three of 
these parameters will be discussed below. The remaining 
two, the initial location and length of the sample, have 
obvious implications for the number density profiles. 



1. Variation in fractional order 7 

The fractional order 7 is a dimensionless quantity 
which defines the degree of the trapping within the 
medium, with a smaller value corresponding to greater 
and longer lasting traps. The maximum value of 7 = 1 
corresponds to "normal transport," which is governed 
by the classical (non-fractional) diffusion advection equa- 
tion. 




0.2 0.4 0.6 O.f 

Normalised Position 
(a) 7 = 1.00 




0.2 0.4 0.6 O.i 

Normalised Position 

(b) 7 = 0.75 




0.2 0.4 0.6 O.i 
Normalised Position 

(c) 7 = 0.50 




0.2 0.4 0.6 0.8 1 

Normalised Position 

(d) 7 = 0.25 

Figure 2: (Color online) Impact of the fractional order 
7 of the trapping distribution on the space-time 
evolution of the number density. In these plots, 
W = 40/L (s-'>') and Dl = (s"^). 
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0.2 0.4 0.6 0.8 1 

Normalised Position 

(a) W = 10.0; Dl = 1.0 




0.2 0.4 0.6 0.8 

Normalised Position 

(a) W = 10.0; Dl = 1.0 




0.2 0.4 0.6 0.8 1 

Normalised Position 

(b) W = 10.0; Dl = 10.0 




0.2 0.4 0.6 0.8 
Normalised Position 

(b) W = 10.0; Dl = 10.0 
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0.2 0.4 0.6 0.8 

Normalised Position 

(c) W = 100.0; Dl = 1.0 




0.2 0.4 0.6 0.8 

Normalised Position 

(c) W = 100.0; Dl = 1.0 
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0.2 0.4 0.6 0.8 1 

Normalised Position 

(d) W = 100.0; Dl = 10.0 

Figure 3: (Color online) Space-time evolution of the 
number density profile for 7 = 0.8. Here, W and Dl are 
normalized to the length of the apparatus and are hence 
both specified in units of s~^. 




0.2 0.4 0.6 0.8 

Normalised Position 
(d) W = 100.0; Dl = 10.0 

Figure 4: (Color online) Space-time evolution of the 
number density profile for 7 = 0.4. Notice that these 
figures use a different time scale to those in Figure [H 
Here, W and Dl are normalized to the length of the 
apparatus and are hence both specified in units of s 
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Figure 5: Experimental time of flight current trace data 
for trinitrofluorenone-polyvinylcarbazole, digitized from 
J5|]. The solid line is the model fit. 



The impact of 7 on the electric current is demonstrated 
in Figure [1] For non-dispersive transport (7 = 1), the 
result is essentially a time independent (displacement) 
current until a sharp cutoff where the charged particles 
exit the system through the electrode. The finite drop 
off time is a reflection of the diffusion in the system. For 
dispersive transport, the departure of the current traces 
from the classical proflles is enhanced as the fractional 
order decreases. The fractional order 7 defines the slopes 
of the two regimes, and hence, characterizes the funda- 
mental shape of the current trace. The relevant relations 
are given in Eq. (jl7|) above. 

Number density profiles corresponding to the afore- 
mentioned current solutions are shown in Figure [2j So- 
lutions for 7 = 1 exhibit a moving Gaussian "pulse" of 
charge carriers, spreading according to Dl and drifting 
according to W. This is shown in Figure 

For 7 < 1, the signature of fractional or dispersive be- 
havior appears. In this mode, the number density profile 
retains a "memory" of the initial sharp spike at z = zq- 
This peak in the density profile does not drift with W, as 
it does in the non-dispersive case. This long persistence 
of the initial condition has previously been mentioned 
in the literature [1, H, [l^. The smaller the value of 7, 
the more dispersive the transport. Indeed, for strongly 
dispersive systems, the spike at z = zq is the most promi- 
nent feature of the entire charge distribution for much of 
its lifetime. This sharp spike is most clearly illustrated 
in the contour plots of Figures [^c] and I2d[ 



2. Impact of the drift velocity W and diffusion coefficient 

Dl 

The fractional drift velocity has units of m/s^ , and 
describes the tendency of the charged particles to drift 
in the positive z direction. The fractional diffusion co- 
efficient has units of m^/s''', and describes the tendency 
of the charged particles to diffuse down the concentra- 
tion gradient. The effects of varying W and Dl are 
demonstrated in Figure [31 for a weakly dispersive sys- 
tem (7 — 0.8); and in Figure U for a strongly dispersive 
system (7 — 0.4). The relevant parameters are indicated 
in the figure captions. For both systems, an increased W 
sweeps the charge carriers further to the right, and an 
increased Dl spreads the swarm over a wider area. 



B. Experimental Results 

To demonstrate the process by which this model may 
be fitted to time-of-fiight experimental data, we consider 
the data for trinitrofiuorenone and polyvinylcarbazole 
(TNF-PVK) presented as Figure 6 of 5]. The data was 
digitized from the scanned plot, and the slopes of the 
two regimes was used to furnish an estimate for 7. We 
used L = 1 to give a normalized length scale; and se- 
lected the initial source location zo to be 0.2, since the 
model is largely insensitive to the location of the source, 
provided it is sufficiently far from the electrodes to avoid 
substantial "back diffusion." 

The intercept of the two straight lines was taken to be 
the transit time ttr, and the following equation was used 
to furnish an estimate of W, which provided a starting 
point for curve fitting: 

the factor of 1/2 being an empirical correction that gives 
better results when compared with the order of magni- 
tude estimate Eq. (|23p. The final remaining parameter 
was initially taken as Dl ~ Vl^/20. 

The parameter estimates discussed above were used as 
the starting point for nonlinear least squares curve fitting. 
The Matlab Curve Fitting Toolbox was used. The result 
of the model fitting is shown in Figure [5] 



V. CONCLUSION 

We have demonstrated a fractional advection diffusion 
equation modeling the hopping transport observed in 
many disordered semiconductors. We have shown that 
the infinite series of Fourier modes [Eq. (|4])] for the 
bounded solution can be collapsed into a closed form 
expression using the Poisson summation theorem [Eq. 
(IT^ ]. It this closed form expression that then facilitates 
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the extraction of model parameters from the experimen- 
tal data using a simple curve fitting routine. We have 
modeled a time of flight experiment by assuming the ini- 
tial condition n{z,to) = noS^z — zo)- We have calculated 
the resultant electric current, and shown that the sum of 
slopes on logarithmic axes is —2, as predicted by other 
models and as verified by experiment. It is possible to 
extend this solution to sources of finite duration or finite 
width, by integrating with respect to to or zq, respec- 
tively. 
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Appendix A: Caputo and Riemann-Liouville forms 
of the Fractional Advection Diffusion equation 

1. Fractional Derivatives 

The two forms of fractional derivative commonly used 
to describe subdiffusive systems are the Caputo deriva- 
tive and the Riemann-Liouville derivative. In what fol- 
lows, we describe fractional partial derivatives with re- 
spect to t in terms of an arbitrary function f{t,x,y, ...). 
For clarify of presentation, the functional dependence of 
/ on the other variables is suppressed, and we write sim- 

ply fit). 

The Caputo derivative of order < a < 1 is defined as 



oD?f{t) ^ 



1 



Til -a) 



(<-r)-"/'Mdr, (Al) 



where /'(r) is the ordinary partial derivative df/dt eval- 
uated at i = T. The Laplace transform of the Caputo 
derivative is: 



e-^' ^Dffit) dt = s'^fis) - s"- V(0), (A2) 



where /(s) is the Laplace transform of fit), and /(O) is 
the initial condition. 

The Riemann-Liouville fractional derivative of order 
< a < 1 is the defined as [201: 



I d 
r(l -a)di Jo 



(i-T)-"/(r)dr. (A3) 



The Laplace transform of a Riemann-Liouville derivative 
is: 



Jo 



where /o is a fractional initial condition: 
f 



fo 



lini / , '^^ \n dT. 



r{l-a)t^oJ^ (t-T) 



(A4) 



2. Fractional Advection-Diffusion Equations 

The first model for dispersive transport was due to 
Scher and Montroll [5] , who used a continuous time ran- 
dom walk (CTRW) where the waiting time probability 
density function has divergent mean. A continuous time 
random walk is characterized by a hopping probability 
density function (pdf) tpiz, t). We consider the decoupled 
case tpiz, t) — A(z)u'(t) where A(z) is the jump length pdf 
and wit) is the waiting time pdf. Under these conditions. 



the CTRW has the Fourier-Laplace space solution |21 



n(fc, s) 



1 — wis) «o(fc) 
~s 1 - A(/c)w(s)' 



(A5) 



where Fourier transformed functions are denoted by ex- 
plicit dependence on the Fourier variable fc, and rioik) is 
the Fourier transformed initial condition. 

We postulate a CTRW where the waiting time pdf has 
divergent mean. Such a pdf has the small s asymptote 
SlH: 



wis) ^ I — (rs)'' 



(A6) 



We further postulate a well-behaved jump length pdf 
with moment generating function 

M^ix)^l + Ahx+^^ + -, 

for first and second moments Mi and M2, respectively. 
This corresponds to a characteristic function (i.e. Fourier 
transform) in the small k limit of: 



A(fc) = Mxiik) - 1 + iMik - 



(A7) 



Substituting these asymptotes into Eq. (jASp . and dis- 
carding terms of order Oiks^) and higher, we obtain: 



nik, s) 



no(fc)s 



,7-1 



St - iWk + DLk^ ' 



(A8) 



where W = Mi/t' and Dl = Mij^T^ . Equation 
(|A8P is the free-space propagator of fractional advection- 
diffusion. By rearranging Eq. \k&\ . one can derive vari- 
ous forms of fractional advection diffusion equation. For 
example, one readily obtains: 



sTn(z, s) - s^-^noiz) + ( W - DLi^ ] n(z, s) = 0, 

(A9) 



— -D — 

dz dz'^ 
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z=0 z=L 




I- z' -I 

L 



Figure 6: (Color online) Simplified time of flight 
schematic used in current derivation. The two 
electrodes at z = and z = L have potentials Vq and 
Vi , respectively. A surface S cuts through the sample at 
z — z'; the volume V is the space between the z = 
electrode and the surface S. 



I=jrj{z',t)dz'+'-^^^{Vo^V^), (Bl) 

where j{z',t) is the conduction current passing through 
the surface S, e is the permittivity of the semiconducting 
material, and A is the area of the electrodes. 

Under typical measuring conditions, the transit time 
ttr is much less than the RC time of the circuit. There- 
fore, we assume that Vq — Vi is essentially constant, and 
then the current is simply the space-averaged conduction 
current: 

I=j £ j{z',t)dz'. (B2) 

The conduction current leaving the volume V is the 
negative rate of change of the charge enclosed: 



which is the Laplace transform of the Caputo fractional 
equation ([T]) . Alternatively, Eq. (|A8p may be rearranged 
to give: 



n[z, sj h s 



which is a fractional integral equation. Inverting the 
Laplace transform in Eq. (jA10[) . and taking an ordi- 
nary partial derivative with respect to time, one obtains 
the following form of the fractional advection diffusion 
equation: 



dn 
'di 



oz oz'^ 



0. 



(All) 



Equation (jAlip is a special case of the fractional 
Fokker-Planck equation [J (2^ , and is equivalent to the 
Caputo fractional advection diffusion equation ([T|) con- 
sidered in this paper. 



Appendix B: Derivation of Current Formula 

Consider a time of flight system where all spatial vari- 
ation is confined to the z direction, normal to the elec- 
trodes. An electrode at z = is held at a potential Vq 
by an external power supply, and the opposite electrode 
a,i z = L has potential Vi and is connected via a resis- 
tor R to the ground, as shown in Figure |6l We define a 
surface S which is normal to the electrodes at a position 
z = z', and a volume V which is the entire area between 
the z = electrode and the surface S. 

The overall current will consist of a conduction current 
and a displacement current. Integrating across the width 
of the sample: 



j{z',t) = -— j qn{z,t)dz. 



Using Eq. p2|) : 



I = — / / n(z,t)dzdz'. 



Changing the order of integration: 



/ = — — — [ [ n(z,t)dz'dz 
L dt Jq 

= ^^^^ {L- z)n{z,t)dz 

= q— I — f zndz — f ndz 
^dt]L.L ./n 



(B3) 



It should be noted that different expressions ex- 
ist within the literature for the current depending on 
whether the paper in question uses a multiple trapping 
model or a hopping model. This is why our current ex- 
pression (jB3p is at first glance not equivalent to the cur- 
rent expressions used by some other authors. Under a 
multiple trapping model, the equivalent is: 



W 



nircc{z,t)dz, 



(B4) 



where nftcc is the distribution of untrapped particles and 
W is the drift velocity of these particles. This formula 
can be obtained by neglecting diffusive fiux to substitute 
j = W^nfroo into Eq. ((B2|) . 
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